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Double beta decay is indispensable to solve the question of the neutrino mass matrix 
together with v oscillation experiments. Recent analysis of the most sensitive experiment 
since nine years - the HEIDELBERG-MOSCOW experiment in Gran-Sasso - yields a 
first indication for the neutrinoless decay mode. This result is the first evidence for lepton 
number violation and proves the neutrino to be a Majorana particle. We give the present 
status of the analysis in this report. It excludes several of the neutrino mass scenarios 
allowed from present neutrino oscillation experiments - only degenerate scenarios and 
those with inverse mass hierarchy survive. This result allows neutrinos to still play an 
important role as dark matter in the Universe. To improve the accuracy of the present 
result, considerably enlarged experiments are required, such as GENIUS. A GENIUS 
Test Facility has been funded and will come into operation by early 2003. 

Keywords: Beta decay, double beta decay; Neutrino mass and mixing; Weak-interaction 
and lepton (including neutrino) aspects. 



1. Introduction 

Double beta decay is the most sensitive probe to test lepton number conservation. 
Further it seems to be the only way to decide about the Dirac or Majorana nature 
of the neutrino. Observation of Ov/3/3 decay would prove that the neutrino is a Ma- 
jorana particle and would be another clear sign of beyond standard model physics. 
Recently atmospheric and solar neutrino oscillation experiments have shown that 
neutrinos are massive. This was the first indication of beyond standard model 
physics. The absolute neutrino mass scale, however, cannot be determined from 
oscillation experiments alone. Double beta decay is indispensable also to solve this 
problem. 

The observable of double beta decay is the effective neutrino mass 
H = | £ U&ni | = \m<£ 1 + \m£) \ + e**> | , 
with U ei denoting elements of the neutrino mixing matrix, m, neutrino mass eigen- 
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states, and 4m relative Majorana CP phases. It can be written in terms of oscillation 
parameters ^ 

|m«| = IC^IV. (!) 



ni 



(2)1 - ITT J 2 . U m 2 



U e2 \'JAm 2 21 +ml (2) 



|C / e 3 lV Am 32 + Am il+ TO ?- ( 3 ) 

The effective mass (m) is related with the half-life for 0v/3f3 decay via 
(-^1/2) ~ ( m;y ) 2 ' ana * f° r ^ ne um it 011 ^1/2 deducible in an experiment we have 



Here a is the isotopical abundance of the /3/3 emitter; M is the active detector 

mass; t is the measuring time; AE is the energy resolution; B is the background 

count rate and e is the efficiency for detecting a /3/3 signal. Determination of the 

17101 

effective mass fixes the absolute scale of the neutrino mass spectrum. IL121. 

In this paper we will discuss the status of double beta decay search. Although 
in the HEIDELBERG-MOSCOW experiment we also have the highest statistics 
(~ 147 000 events) for 2z//3/3 decay (see EI and ®) we shall concentrate in this paper 
on the ncutrinolcss decay mode. We shall, in section 2, discuss the recent evidence 
for the neutrinoless decay mode, from the HEIDELBERG-MOSCOW experiment, 
and the consequences for the neutrino mass scenarios which could be realized in 
nature. In section 3 we discuss the possible future potential of Qvf3(3 experiments, 
which could improve the present accuracy. 



2. Evidence for the Neutrinoless Decay Mode 

The status of present double beta experiments is shown in Fig. ^ and is exten- 
sively discussed in 1*1 The HEIDELBERG-MOSCOW experiment using the largest 
source strength of 11 kg of enriched 76 Ge (enrichment 86%) in form of five HP Ge- 
detectors is r unning since August 1990 in the Gran-Sasso underground laboratory 
M4 | 5JlKj6J32 | anc j j g smce nme y ears now the most sensitive double beta experiment 
worldwide. 



2.1. Data from the HEIDELBERG-MOSCOW Experiment 

The data taken in the period August 1990 - May 2000 (54.9813 kg y, or 723.4 4 mol- 
years) are shown in Fig-Oin the section around the Qpp value of 2039.006 
Fig. |21 is identical with Fig. ^ in ^ except that we show here t he or iginal energy 
binning of the data of 0.36 keV. These data have been analysed l ^ l ^ l ^ l with various 
statistical methods, with the Maximum Likelihood Method and in particular also 
with the Bayesian method (see, e.g. This method is particularly suited for 
low counting rates, where the data follow a Poisson distribution, that cannot be 
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Fig. 1. Present sensitivity, and expectation for the future, of the most promising (3/3 experiments. 
Given are limits for (m), except for the HEIDELBERG-MOSCOW experiment where the recently 
observed value is given (95% c.l. range and best value). Framed parts of the bars: present status; 
not framed parts: future expectation for running experiments; solid and dashed lines: experiments 
under construction or proposed, respectively. For references see l ^ l ^ | 4 | 39 | 50 | 
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Fig. 2. The spectrum taken with the 76 Ge detectors Nr. 1,2,3,4,5 over the period August 1990 
- May 2000 (54.9813 kg y) in the original 0.36 keV binning, in the energy range 2000 - 2100 keV. 
Simultaneous fit of the 214 Bi lines and the two high-energy lines yield a probability for a line at 
2039.0 keV of 91%. 
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appro ximated by a Gaussian. Details and the results of the analysis are given in 
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Fig. 3. Top: Probability K that a line exists at a given energy in the range of 2000-2080 keV 
derived via Bayesian inference from the spectrum shown in Fig. [2] Bottom: Result of a Bayesian 
scan for lines as in the left part of this figure, but in an energy range of ±5cr around 1 



Our peak search procedure (for details see EG}) reproduces (see ^^ j) 7-lines at 
the positions of known weak lines from the decay of 214 Bi at 2010.7, 2016.7, 2021.8 
and 2052.9 keV^l In addition, a lin e cente red at 2039 keV shows up (see Fig. EJ. 
This is compatible with the Q- value l ^ l ^ l of the double beta decay process. The 
Bayesian analysis yields, when analysing a ±5er range around Qpp (which is the 
usual procedure when searching for resonances in high-energy physics) a confidence 
level (i.e. the probability K) for a line to exist at 2039.0 keV of 96.5 % c.l. (2.1 
a) (see Fig.OJl. We repeated the analysis for the same data, but except detector 4, 
which had no muon shield and a slightly worse energy resolution (46.502 kg y). The 
probability we find for a line at 2039.0 keV in this case is 97.4% (2.2 a) l 1 ! ^ ! 4 ! 

Fitting a wide range of the spectrum yields a line at 2039 keV at 91% c.l. (see 
Fig©. 

We also applied the Feldman-Cousins method 031 This method (which does 
not use the information that the line is Gaussian) finds a line at 2039 keV on a 
confidence level of 3.1 a (99.8% c.l.). In addition to the line at 2039 keV we find 
candidates for lines at energies beyond 2060 keV and around 2030 keV, which at 
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present cannot be attributed. This is a task of nuclear spectroscopy. 

Important further information can be obtained from the time structures of the 
individual events. Double beta events should behave as single site events (see Fig. 
El top), i.e. clearly different from a multiple scattered 7-event (see Fig. El bottom). 
It is possible to differentiate between these different types of event s by pulse shape 
analysis. We have developped three methods of pulse shape analysis ^ W during 
the last seven years, one of which has been patented and therefore only published 
recently. 
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Fig. 4. Top: Shape of one candidate for 0vf3/3 decay classified as SSE by all three methods 
of pulse shape discrimination. Bottom: Shape of one candidate classified as MSE by all three 
methods. 

Installation of Pulse Shape Analysis (PSA) has been performed in 1995 for the 
four large detectors. Detector Nr. 5 runs since February 1995, detectors 2,3,4 since 
November 1995 with PSA. The measuring time with PSA from November 1995 
until May 2000 is 36.532 kg years, for detectors 2,3,5 it is 28.053 kg y. 

In the SSE spectrum obtained under the restriction that the signal simultane- 
ously fulfills the criteria of all three methods for a single site event, we find again 
indication of a line at 2039. OkeV (see Figs.EUHJ. 





February 1, 2008 20:45 WSPC/INSTRUCTION FILE HVKK-StonBr02- 
corr 



6 H. V. KLAPDOR-KLEINGROTHA US 




energy [keV] 



Fig. 5. Scan for lines in the single site event spectrum taken from 1995-2000 with detectors Nr. 
2,3,5, (Fig. 0, with the Bayesian method. Top: Energy range 2000 -2080 keV. Bottom: Energy 
range of analysis ±4.4cr around Q/3/3. 
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Fig. 6. Sum spectrum of single site events, measured with the detectors Nr. 2,3,5 operated with 
pulse shape analysis in the period November 1995 to May 2000 (28.053 kgy), summed to 1 keV 
bins. Only events identified as single site events (SSE) by all three pulse shape analysis methods 
have been accepted. The curve results from Bayesian inference in the way explained in 
sec. 3. When corrected for the efficiency of SSE identification (see text), this leads to the following 
value for the half-life: T°"; 2 =(0.88 - 22.38) X 10 25 y (90% c.l.). 
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We find 9 SSE events in the region 2034.1 - 2044.9keV (± 3cr around Q^). 
Bayes analysis of the range 2032 - 2046 keV yields a signal of single site events, as 
expected for neutrinoless double beta decay, with 96.8% c.l. at the Qpp value. The 
Feldman-Cousins method gives a signal at 2039.0 keV of 2.8 a (99.4%). 

The analysis of the line at 2039.0 keV before correction for the efficiency yields 
4.6events (best value) or (0.3 - 8.0) events within 95% c.l. ((2.1 - 6.8) events within 
68.3% c.l.). Corrected for the efficiency to identify an SSE signal by successive appli- 
cation of all three PSA methods, which is 0.55 ± 0.10, we obtain a Oi//3/3 signal with 
92.4% c.l.. The signal is (3.6 - 12.5) events with 68.3% c.l. (best value 8.3events). 
Thus, with proper normalization concerning the running times (kgy) of the full 
and the SSE spectra, we see that almost the full signal remains after the single site 
cut (best value), while the 214 Bi lines (best values) are considerably reduced. We 
have used a 238 Th source to test the PSA method. We find the reduction of the 
2103 keV and 2614 keV 228 Th lines (known to be multiple site or mainly multiple 
site), relative to the 1592 keV 22S Th line (known to be single site), shown in Fig.0 
This proves that the PSA method works efficiently. Essentially the same reduction 
as for the Th lines at 2103 and 2614 keV and for the weak Bi lines is found for the 
strong 214 Bi lines (e.g. at 609.6 and 1763.9 keV (Fig. 01). 

The possibility, that the single site signal is the double escape line corresponding 
to a (much more intense!) full energy peak of a 7-line, at 2039+1022=3061 keV is 
excluded from the high-energy part of our spectrum. 

' * Range of Bi lines + Ovpp Signal Calibration Lines 




2010 2016 2021 2053 2039 609.6 1763.9 1592 2103 2614 

' v / ' ' " ' 



Fig. 7. Relative suppression ratios: Remaining intensity after pulse shape analysis compared to 
the intensity in the full spectrum. Right: Result of a calibration measurement with a Th source - 
ratio of the intensities of the 1592 keV line (double escape peak, known to be 100% SSE), set to 
1. The intensities of the 2203 keV line (single escape peak, known to be 100% MSE) are strongly 
reduced (error bars are ±lo\ The same order of reduction is found for the strong Bi lines occuring 
in our spectrum - shown in this figure are the lines at 609.4 and 1763.9 keV. Left: The lines in the 
range of weak statistics around the line at 2039 keV (shown are ratios of best fit values). The Bi 
lines are reduced compared to the line at 2039 keV (set to 1), as to the 1592 keV SSE Th line. 
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A very careful simulation of the different components of radioactive background 
in the Heidelberg-Moscow experim ent h as been performed recently by a new Monte 
Carlo program basing on GEANT4EEDI. This simulation uses a new event generator 
for simulation of radioactive decays basing on ENSDF-data and describes the decay 
of arbitrary radioactive isotopes including a, [3 and 7 emission as well as conversion 
electrons and X ray emission. Also included in the simulation is the influence of 
neutrons in the energy range from thermal to high energies up to 100 MeV on the 
measured spectrum. Elastic and inelastic reactions, and capture have been consid- 
ered, and the corresponding production of radioactive isotopes in the materials of 
the setup. The neutron fluxes and energy distributions were taken from published 
measurements performed in the Gran Sasso. Also simulated was the influence of 
the cosmic muon flux measured in the Gran Sasso, on the measured spectrum. 

The simulation gives no indication that the signal at 2039 keV comes from a 
known background line. In particular, the simulation shows, that e.g. decays of 
77 Ge, 76 Ga or 228 Ac, should not lead to signals visible in our measured spectra 
near the signal at Qpp. For details we refer to^l 



2.2. Comparison with earlier results 

We applied the same methods of peak search as used i n ou r analysis to the spectrum, 
measured in the Ge experiment by Caldwell et al. ^2 more than a decade ago. 
These authors had the most sensitive experiment using natural Ge detectors (7.8% 
abundance of 76 Ge). With their background being a factor of 9 higher than in the 
present experiment, and their measuring time of 22.6 kg y, they have a statistics for 
the background larger by a factor of almost 4 in their (very similar) experiment. 
This allows helpful conclusions about the nature of the background. 

The peak scanning finds El (Fig. [HJ) indications for peaks essentially at the same 
energies as in Fig. EI This shows that these peaks are not fluctuations. In particular 
it sees the 2010.78, 2016.7, 2021.6 and 2052.94keV 2U Bi lines. It finds, however, 
no line at Qpp. This is consistent with the expectation from the rate found from 
the HEIDELBERG-MOSCOW experiment. About 16 observed events in the latter 
correspond to 0.6 expected events in the Caldwell experiment, because of the use of 
non-enriched material and the shorter measuring time. 

The first experiment using enriched ( but no t high-purity) 76 Ge detectors per- 
formed by Kirpichnikov and coworkers l ^ l ^ l because of their low statistics of 
2.95 kg y would expect 0.9 counts. Their result is consistent with this expectation. 
Another Ge experiment (IGEX) using 8.8 kg of enriched 76 Ge, but collecting since 
beginning of the ex peri ment in the early nineties till shutdown in end of 1999 only 
8.8 kg y of statistics 1*2^ could expect, according to our result, about 2.6 events. The 
result of that measurement is also consistent with the expectation. 
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Fig. 8. Peak scanning of the spectrum measured by Caldwell et al. iJd^ with the Maximum 
Likelihood method (upper part), and with the Bayesian method (lower part) (as in Figs. I3l6l 
(selO). 



2.3. Proofs and Disproofs 

The result described in section 2.1 has been questioned in some papers (Aalseth et 
al, |hep-ex/0202018| Feruglio et al., Nucl. Phys. B 637(2002)345; Zdesenko et al, 
Phys. Lett. B 546(2002) 206). We think that we have shown in a convincing way 
that these claims against our results are incorrect in various ways. In particular 
the estimates of the intensities of the 214 Bi lines in the first two papers do not 
take into account the effect of true coincidence summing, which can lead to drastic 
underestimation of the intensities. A correct estimate would also require a Monte 
Carlo simulation of our setup, which has not been performed in the above papers. 
All of these papers, when discussing the choice of the width of the searc h w indow, 
seem to ignore th e resu lts of the statistical simulations we published inEEE. For 
details we refer to^EEl 



2.4. Half-Life and Effective Neutrino Mass 

Having shown that the signal at Qpp consists of single site events and is not a 
7-line, we translate the observed number of events into half-lifes. We obtain ID^J4J 
Tfl 2 = (0.8 - 18.3) x 10 25 y (95% c.l.) with a best value of 1.5 x 10 25 y Assuming 
that the 0^/3/3 amplitude is dominated by the neutrino mass mechanism, we obtain, 
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with the nuclear matrix element froml^an effective mass of (m !/ )=(0.11 - 0.56) eV 
(95% c.l.). 

The result obtained is consistent with all other double beta experiments - 
which still reach less sensitivity. The most sensitive experiments following the 
HEIDELBERG-MOSCOW experiment are the geochemical 128 Te experiment with 
T° v /2 > 2(7.7) x 10 24 y (68% c.l.),Elthe 136 Ae experiment by the DAMA group with 

T°"> 1.2 x 10 24 y (90% c.l.jSO a second 76 Ge experiment with T°J 2 > 1.2 x 10 24 
y and a nat Ge experiment with T"j 2 > 1 x 10 24 y EH21 Other experiments 
are already about a factor of 100 less sensitive concerning the 0v/3f3 half-life: the 
Gotthard TPC experiment with 136 Xe yields^ T°v > 4.4 x 10 23 y (90% c.l.) and 

the Milano Mibeta cryodetector experiment 4 8 T°J 2 > 1.44 x 10 23 y (90% c.l.). 

Another experiment ^Bl with enriched 76 Ge, which has stopped operation in 
1999 after reaching a significance of 8.8 kg y, yields (if one believes their method of 
'visual inspection' in their data analysis), in a conservative analysis, a limit of about 
T^ 2 > 5 x 10 24 y (90% c.l.). The 128 Te geochemical experiment yields (m„) < 1.1 

eV (68 % c.l.) EH, the DAMA 136 Ae experiment (m„) < (1.1 - 2.9)eV^and the 
130 Te cryogenic experiment yields {m u ) < 1.8 eV l^al 

Concluding we obtain, with about 95% probability, first evidence for the neutri- 
noless double beta decay mode. As a consequence, at this confidence level, lepton 
number is not conserved. Further the neutrino is a Majorana particle. If the Qv(3(3 
amplitude is dominated by exchange of a massive neutrino the effective mass (m) 
is deduced to be (m) = (0.11 - 0.56) eV (95% c.l.), with best value of 0.39eV. Al- 
lowing conservatively for an uncertainty of the nuclear matrix elements of ±50% 
(for detailed discussions of the status of nuclear matrix elements we refer to 
and references therein) this range may widen to (m) = (0.05 - 0.84) eV (95% c.l.). 

Assuming other mechanisms to dominate the Qv(5f3 decay amplitude, the re- 
sult allows to set stringent limits on parameters of SUSY models, leptoquarks, 
compositeness, masses of heavy neutrinos, the right-handed W boson and possible 
violation of Lorentz invariance and equivalence principle in the neutrino sector. For 
a discussion and for references we refer to W I ^7 l sM I 38 l v>Q I 

With the limit deduced for the effective neutrino mass, the HEIDELBERG- 
MOSCOW experiment excludes several of the neutrino mass scenarios allowed from 
present neutrino oscillation experiments (see Fig. EJ) - allowing only for degenerate 
and inverse hierarchy mass scenarios^. 

Assuming the degenerate scenarios to be realized in nature we fix - according to 
the formulae derived in^l- the common mass eigenvalue of the degenerate neutrinos 
to m = (0.05 - 3.4) eV. Part of the upper range is already excluded by tritium 
experiments, which give a limit of m < 2.2-2.8 eV (95% c.l.)' 4 ^'. The full range can 
only partly (down to ~ 0.5 eV) be checked by future tritium decay experiments, 
but could be checked by some future (3/3 experiments (see, e.g., next section). The 
deduced best value for the mass is consistent with expectations from experimental 
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HEIDELBERG- MOSCOW, FIRST OvPP Evidence 
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Fig. 9. The impact of the evidence obtained for neutrinoless double beta decay in this paper 
(best value of the effective neutrino mass (m) = 0.39 eV, 95% confidence range (0.05 - 0.84) eV - 
allowing already for an uncertainty of the nuclear matrix element of a factor of ± 50%) on possible 
neutrino mass schemes. The bars denote allowed ranges of (m) in different neutrino mass scenarios, 
still allowed by neutrino oscillation experiments (seelHI). Hierarchical models are excluded by the 
new 0vf3f3 decay result. Also shown are the expected sensitivities for the future potential double 
beta experiments CUORE, MOON, EXO and the 1 ton and 10 ton project of GENIUSES 38 ! 51 ! 29 ! 



fj, — ► cy branching limits in models assuming the generating mechanism for the 
neutrino mass to be also responsible for the recent indication for as anomalous 
magnetic moment of the muon ESI. It lies in a range of interest also for Z-burst 
models recently discusse d as expl anation for super-high energy cosmic ray events 
beyond the GKZ-cutoff A recent model with underlying A4 symmetry 

for the neutrino mixing matrix also leads to degenerate neutrino masses consistent 
with the present result from Oz//3/3 decay ESI xhe range of (to) fixed in this work 
is, alread y now , in the range to be explored by the satellite experiments MAP and 
PLANCK EB3 (see Fig. EH . 

The neutrino mass deduced leads to 0.002> Vl u h? < 0.1 and thus may allow 
neutrinos to still play an important role as hot dark matter in the Universe E2 

3. Future of (3(3 Experiments - GENIUS and Other Proposals 

With the HEIDELBERG-MOSCOW experiment, the era of the small smart experi- 
ments is over. New approaches and considerably enlarged experiments (as discussed, 
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mi (eV) HEIDELBERG-MOSCOW Positive EVIDENCE 




Fig. 10. Double beta decay observable (m) and oscillation parameters: The case for degenerate 
neutrinos. Plotted on the axes are the overall scale of neutrino masses mo and mixing t an 2 9 19 . 
Also shown is a cosmological bound deduced from a fit of CMB and large scale structure and 
the expected sensitivity o f th e satellite experiments MAP and PLANCK. The present limit from 
tritium f} decay of 2.2 eV 

34 

would lie near the top of the figure. The range of (m) fixed by the 
HEIDELBERG-MOSCOW experiment d is, in th e ca se of small solar neutrino mixing, already in 
the range to be explored by MAP and PLANCK EH 



e.g. [ n i a2iQUlsUl^li^s2^ w iH be required in future to fix the neutrino mass with higher 
accuracy. 

Since it was realized in the HEIDELBERG-MOSCOW experiment, that the re- 
maining small background is coming from the material close to the detector (holder, 
copper cap, ...), elimination of any material close to the detector will be decisive. 
Experiments which do not take this into account, like, e.g. CUORE^I and MAJO- 
RANA ESI will allow at best only rather limited steps in sensitivity. Furthermore 
there is the problem in cryodetectors that they cannot differentiate between a (3 
and a 7 signal, as this is possible in Ge experiments. 

Another crucial point is - see eq. (4) - the energy resolution, which can be 
optimized only in experiments using Germanium detectors or bolometers. It will be 
difficult to probe evidence for this rare decay mode in experiments, which have to 
work - as result of their limited resolution - with e nerg y windo ws a round of 
several hundreds of keV, such as NEMO IlP3 Exd^CAMEcPl. 

Another important point is (see eq. 4), the efficiency of a detector for detection 
of a (3(3 signal. For example, with 14% efficiency a potential future 100 kg 82 Se 
NEMO experiment would be, because of its low efficiency, equivalent only to a 
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10kg experiment (not talking about the energy resolution). 

In the first prop osal for a third generation double beta experiment, the GE- 
NIUS proposal 1 ^7 1 ^8 1 3^ 1 5 1 1 ^9 1 the idea is to use 'naked' Germanium detectors in a 
huge tank of liquid nitrogen. It seems to be at present the only proposal, which can 
fulfill both requirements mentioned above - to increase the detector mass and si- 
multaneously reduce the background drastically. GENIUS would - with only 100 kg 
of enriched 76 Ge - increase the confidence level of the present pulse shape dis- 
criminated Qv(3(3 signal to 4er within one year, and to 7a within three years of 
measurement (a confirmation on a 4cr level by the MAJORANA project would at 
least need ^230 years, the CUORE project would need (ignoring for the moment 
the problem of identification of the signal as a (3(3 signal) 3700 years). With ten 
tons of enriched 76 Ge GENIUS should be capable to investigate also whether the 
neutrino mass mechanism or another mechanism (see, e.g. ™) is dominating the 
Qu(3(3 decay amplitude. A GENIUS Test Facili ty is a t present under construction 
in the GRAN SASSO Underground Laboratory BM3 

4. Conclusion 

The status of present double beta decay search has been discussed, and recent evi- 
dence for a non- vanishing Majorana neutrino mass obtained by the HEIDELBERG- 
MOSCOW experiment has been presented. The latter opens a new era in space-time 
structure 1 ^. It has been shown that the Majorana nature of the neutrino tells 
us that spacetime does realize a construct that is central to construction of super- 
symmetric theories. 

Future projects to improve the present accuracy of the effective neutrino mass 
have been briefly discussed. The most sensitive of them and perhaps at the same 
time most realistic one, is the GENIUS project. GENIUS is the only of the new 
projects which simultaneously has a huge potential for cold dark matter s earch, and 
for real-time detection of low-energy neutrinos (see l ^7 | 33 | 35 | 38|4 2 39 o0j 
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